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Abstract 

The quantum hadrodynamics are studied by using non-perturbative renormal¬ 
ization group equations. Approximate discrete 75 -symmetry is studied. It is 
shown that the contributions of one-loop diagrams are important for effective 
potential. A relation between the local-potential approximation and the usual 
Hartree approximation is discussed. The local-potential approximation includes 
contribution of the Fock term to nucleon self-energy, contribution in the usual 
random phase approximation to meson self-energy and vertex corrections in part 
as well as the Hartree contribution to the effective potential. Approximate dis¬ 
crete 75 -symmetry is studied phenomenologically in the Hartree approximation. 
This approximate symmetry may exist in nuclear matter. 
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1 Introduction 


In the past two decades, nuclei and nuclear matter have been studied in the 
framework of quantum hadrodynamics (QHD).^, ^ The meson mean-field the¬ 
ory of nuclear matter Q] has produced successful results to account for the satu¬ 
ration properties at normal nuclear density. Following these successes, many 
studies and modifications have been performed on relativistic nuclear mod¬ 
els. One of these modifications is the inclusion of vacuum fluctuation effects, 
which cause divergences in physical quantities when they are naively calculated. 
Chin[^ estimated vacuum fluctuation effects in the Hartree approximation by 
using a renormalization procedure, and found that vacuum fluctuation effects 
make the incompressibility of nuclear matter smaller and closer to the empirical 
value than in the original Walecka model. 

Although the relation between QHD and the underlying fundamental theory, 
i.e., QCD, is an open question, it is natural that QHD is not valid at very high- 
energy. In this point of view, a cutoff or a form factor should be introduced into 
the theory of QHD. One may introduce the cutoff[|j or the form factor ||] to avoid 
the instability of the meson propagators in the random phase approximation 
(RPA)I^. Cohen|Q introduced a four-dimensional cutoff into the relativistic 
Hartree calculation and found that the vacuum energy contribution may be 
somewhat different from that in the ordinary renormalization procedures, if the 
cutoff is not so large. 

One troublesome problem in use of the cutoff is that physical results depend 
on a value of the cutoff and a shape of the regulator which are introduced 
into the theory by hand. It is difficult to determine a suitable value of cutoff 
and a appropriate shape of the regulator phenomenologically. For the high- 
energy physics near the limitation of the theory of the particle physics, Lepage^ 
proposed to use the extended effective action (Lagrangian) to eliminate such a 
dependence on cutoff which is by hand introduced into the theory, rather than 
to search the phenomenol^ically favorable cutoff and regulator (or form factor). 
In the series of papers) , we have studied nuclear matter properties and 

the vertex corrections in the framework of the cutoff field theory, eliminating 
the cutoff dependence by using the Lepage’s method. 

The extended action in that method is essentially the same as the one in the 
non-perturbative renormalization group (NPRG) equation. [0 0 0, 0, 0, 
0000, 0, 0 In fact, the former is the perturbative approximation for 
the later. Therefore, it is natural to use the NPRG equations in QHD (more 
generally in nuclear physics), since the QHD has strong couplings. 

In this paper, we study QHD using the NPRG equation. For simplicity, 
we restrict our discussion to the one at the zero-density. We also only treat 
the nucleon and the cr-meson in the NPRG equation. Extension to the finite 
baryon density and inclusion of the other meson are the problems in the furture. 
This paper is organized as follows. In §2, we derive the NPRG equations for 
the tj-nucleon system. In §3, we study the discrete 75 -symmetry by solving the 
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NPRG equations in the local-potential approximation. 0 , |il> 0 S 0 

It is shown that approximate 75 -symmetry may be preserved to some extent 
in analogy of the exact symmetry. 118 In §4, the relation between the local- 
potential approximation and the traditional Hartree approximation in QHD is 
discussed. The local-potential approximation includes contribution of the Fock 
term to nucleon self-energy, contribution in the random phase approximation 
(RPA) to meson self-energy and vertex corrections in part as well as the Hartree 
contribution to the effective potential. In §5, the approximate 75 -symmetry is 
studied phenomenologically by using the Hartree approximation. It is shown 
that this approximate symmetry may exist in nuclear matter. Section 6 is 
devoted to a summary. 


2 Non-perturbative renormalization group equa¬ 
tions 

in quantum hadrodynamics 


We use the NPRG equations with a sharp cutoff function 0(A — 

^ 1^, ^ 1^ where A is cutoff and p is four momentum in Euclidean space. 
Using this cutoff function, the Euclidean action ^[G, A] is given by 


S[e;A| = y:l 


d^pi 


d*Pn 


{2tt)‘^S‘^{pi + ... +Pn) 


„ n\ J (27r)4'"7 (27r)4 

x6»(A - |pi|)...6»(A - \pn\)g^^...i„lJ>l...pn;A)Qi^i-Pl)■..<^^^{-Pn), (1) 


where 0i G {(/), V'li tp are the nucleon field and cr-meson field and g de¬ 

notes a coupling, respectively. The generating functional Z[J] for the Schwinger 
functions (the Euclidean-space Green functions) is given by 


Z[J] = j [d0], 


exp (—S'[0; A]-I- / d^xJi{x)Qi{x) 


( 2 ) 


where Ji are sources to the field 0 ^. 

Next we put A = Aq exp (—t) where Aq is the initial value of the cutoff and t 
is a dimensionless parameter which runs from 0 to infinity. The NPRG equations 
are derived by varying the cutoff and by considering the response of the theory 
subject to the constraint that the Schwinger functions remain unaltered. (For 
the details of the derivation, see Ref. 18). ) It is given by a differential equation 
with respect to t as follows. 



1 1 

2 ^ 

1 1 

2jt 


' 1 — 


dPp 


str lnSij(p, -p) 


'l-5t<|35|,|g|<l 


d^p d^q ,^, 75S[0;t] ^_i SS[Q-,t] 
( 2^)4 ( 2 ^) 4 ^ ^ SQ^^p) *4 ^P’^>SQ^{q) 
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+ 

X 


4-E 

i 


ipi<i 


d'^p 


0i(-p) de 


iBi + p^ 


dpf^ 


5 

SQi{p) 


( 3 ) 


where and jq. are the naive and the anomalous dimension of field 0i(a;), 
respectively, and hi is zero if 0 ^ is a boson field and unity if 0 ^ is a fermion 
held. The symbol ”str” in equation (i denotes supertrace. We remark that 
the four momentum and the held in Eq. § have been rescaled and become 
dimensionless by the replacement 

^ T’ ^ ^ei(A)0*(p), (4) 

where Dq^ (A) is given by following equation 

AEini:)e.(A) = 4-de,+70^. (5) 

Correspondingly, the couplings have been rescaled and become dimensionless 
by the replacement 

n 

5*i...z„(pi,...,Pn;A) = A^""‘‘5,i...,„(pi,...,p„;A)]^i:)Q^ (A). (6) 

A A ij 

i=i 


Using these dimensionless quantities, S{Q;t) is dehned by 

S{Q;t) = S{Q;A) 

The inverse propagators Sij (p, q) are given by 


S^j{p, q) 


5QMdQ^{py 


(7) 

( 8 ) 


where Hij is unity if both the held 0^ is an even element of the Grassmann al¬ 
gebra and the held Oj is an odd element of the Grassmann algebra and vanishes 
otherwise. 

In this section, below(except for (14)), we work only with dimensionless 
helds and momentum vectors and drop the caret notation on all quantities for 
the simplicity of the notation. 

Equation i) is exact for all quantum corrections including any higher or¬ 
der correction in perturbation theory. However, it is the functional differential 
equation and difficult to be solved. To solve it, we ignore all interactions involv¬ 
ing derivative couplings and furthermore set the wavefunction renormalizations 
to unity, therefore, the anomaous dimension 70 ^ is neglected. This is called 
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the local-potential (or action ) approximation 0 , 0 0 0 , 0 , 0 and the 
lowest-term of a systematic momentum (or derivative) expansion. The locally 
approximated action is given by 

S[e = const.; A] = {2Tr)^6^{0)U{Q; t), (9) 


where a generalized potential is given by 


[/[0;t]=^a,....„(t)0.,...0.„, (10) 

n 

Qi is defined by 

0,(p) = {2T,f5\p)Q,. (11) 

and denotes the coupling which does not depend on p. Below we do 

not denote the t-dependence of the couplings explicitly for the simplicity 

of the notation. 

In this approximation, Eq. (^ is reduced to 


d , 1 1 


'l-i5t<|p|<l 


strln[r.,(p) + (-)'‘---^[/(0,t)] 


(27r) 


90^90, 
d 


+ [4-^de.02^]C/(0,i)- (12) 


r,,(p) in Eq. 0) is a kinetic term and given by 


T{p) 


0 

0 


0 \ 

0 - 17 ^ • p 0 /i , 

-iy - p^ h 0 I 


(13) 


where 7 is the 7 -matrix in Euclidean space, h is the unit matrix in isospin 
space, and the symbol ® denotes a direct product. We remark that, in the 
local-potential approximation, the second replacement in Eqs. (^) and the re¬ 
placement (^ are reduced to 

0, = A-''«.0, and (14) 

n 

where = 4 — is a naive dimension of the coupling without 

Z = 1 

the caret. In Eq. (|0), we see that the r. h. s. has a one-loop structure 
with ’’full propagators” which include the quantum corrections in higher energy 
region. 

Furthermore, if we assume the generalized potential {7(0, t) depends only 
on the real scalar field (j) and a = ipiil’i (we remark that the sum is taken over 
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both the spinor and the isospin indices in a), Eq. ( 0 ) is reduced to 

a;t) - iaU^{4>,a\t) 

+ Y^ln(l + - ^ln(l + C/J((/), cr;t)) 


( 15 ) 


where 


11 = 2 - 


+ C/.0( 


1 


1 + C^00 




(16) 


A is the degree of freedom of isospin of fermion (in nucleon system, A = 2) and 
C/cr = -^U etc.. As is seen in §3, the first, the second and the third terms in the 
right-hand side (r. h. s.) of Eq. ( p^ reflects the effects of the naive dimensions 
of the original couplings without the caret. Below, we call this part ’’naive 
dimensional part”. The fourth and the fifth terms in the r. h. s. of Eq. 
look like a bosonic and a fermionic one-loop contributions, respectively. Below, 
we call them ’’bosonic part” and ’’fermionic part”, respectively. The last term 
in r. h. s. of Eq. di) has both bosonic and fermionic contributions. We call 
this part ’’boson-fermion (BE) part”. As is seen in (|l5|), the BE part has at 
least one cr. Therefore, this term contributes directly only to the evolutions of 
the couplings of the the interactions which have at least one a. 

Now, the functional differential equation (^) is reduced to the partial dif¬ 
ferential equation (jl^). In the next section, we study approximate discrete 
75 -symmetry using Eq. (15). 


3 Approximate discrete 75 -symmetry in quan¬ 
tum hadrodynamics 

In this section, we use the coupling constant Ci^m which is defined by 

00 

?/(</), (T;t) = ^ a™ (17) 

0 

rather than ones defined in Eq. (0). In the actual calculation, we truncate 
infinite series in Eq. as follows. 


8 5 2 

U{^, a-t)=Yl ^ ^ (18) 

;=o 1=0 1=0 

This means that we neglect the higher order interactions whose naive 

dimensions are greater than 8 . If we put (^^ into the first three terms of the r. 
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h. s. of Eq. dil), we get 

8 5 2 

^(4 - + ^(1 - l)Ci,,(l)‘a + ^(-2 - /)Q,2</'V2. 

1^0 1^0 1^0 

8 5 2 

= ^ dl,Ci,ocl^‘ + Y, dl,Ci,ic^^a + Y dUCi,2<t>'<y^- (19) 

1=0 1=0 1=0 

We see that the naive dimension df of the couplings Ci^m appears before the 
corresponding coupling. 

In high energy physics, a discrete 75 -symmetry is sometimes imposed on the 
generalized potential, [Q i.e, U{(j),a;t) is invariant under the transformation, 

</>' = -<!), V'' = = TV’Ts, cr' = -a, ( 20 ) 


where the plus (minus) sign in the second equation and the minus (plus) sign 
in the third equation are for proton (neutron). If we extend the theory to 
include pion, this discrete symmetry may be embedded in the continuous chiral 
symmetry. The discrete 75 -symmetry forbides the ’’odd” interaction terms with 
l + m equal an odd integer. Therefore, explicit fermion mass terms are forbidden 
as in the chiral symmetric theory. 

However, in QHD, the nucleon already has a explicit mass term, so we do 
not impose this symmetry on U{4),a\t). However, if we impose this symmetry 
on the initial conditions for the differential equation, i.e., put all of the ’’odd” 
coupling Ci^rn with l + m equal an odd integer to zero at t = 0 , they are remains 
zero through the whole evolution. This results is understood as follow. If we 
expand the Eq. ( 0 ) in power of cj) and cr, each term of the expansion has a form 
as 


l 


(j) a 


dCl,r 

dt 






{l + 2C2,oy^il + Cl,)>^^^ 


( 21 ) 


where ji and ki are nonnegative integers, and Ai is a constant which is indepen¬ 
dent of the fields and couplings. The sum of the r. h. s. in Eq. ( 0 ) correspond 
to the variation of the coupling Ci^m in the infinitesimal interval St. It is easy 
to observe that the sum 

Tli 

'^{Isi+rus,) ( 22 ) 

is an odd integer if l+m is an odd integer and is an even integer if l+m is an even 
integer. This means that, in the case of the odd coupling evolutions equation, 
the numerators in the sum of Eq. (|^ ) have an odd number of odd couplings. 
Therefore, the numerators have at least one odd coupling. For example, a (j)^- 
term in the expansion of the bosonic part 


1 

IGtt^ 


ln(l -I- 


(23) 
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is 


1 


(3!)(167r2) 


120 


Q 


5,0 


1 + 2 ( 72,0 


-432 


C 3 ,0^74^0 

(1 + 2 C 2 , o )2 


432 


^3,0 


(1 + 2 ( 72,0 



(24) 


The numerators in each term of Eq. (|24| ) have at least one odd coupling, there¬ 
fore Eq. is zero at t = 0 if all odd couplings are zero at t = 0. The other 
coefficients of in the r. h. s. of (H also become zero at t = 0. Therefore, 
we get 


dC3,o 

dt 


= 0 


(25) 


at t = 0. Samely, the other odd couplings are also unchanged. Therefore, 
the odd couplings remain zero through the whole evolution if all of the odd 
couplings are put to be zero at t = 0. This mechanism may raise the possibility 
that approximate discrete 75 -symmetry remains through the evolution as well 
as the exact symmetry. Below, we examine this possibility. 

In hgs. we show some results with the initial conditions with ( 72,0 = 

0.05, ( 74,0 = 2.5, (7o,i = 0.3, (7i,i = —2.0 and the other couplings equal to zero 
at t = 0. ( We remark that we use the negative value of (7i,i to correspond to 
the usual definition of Yukawa coupling in QHD.|Q See also Eq. ji^ ) in §4.) 
In this case, only the fermion mass term (7o,icr break the discrete 75 -symmetry 
at the initial. Although, in this paper, we show only the results with one set 
of initial conditions, a qualitative features are similar if we change the initial 
conditions to some extent. 

It is seen that the couplings with the positive naive dimension rapidly in¬ 
crease in the large t region, while the coupling with the negative naive dimension 
rapidly decrease in the same region of t. The couplings whose naive dimensions 
are zero seem to approach a constant value in the large t limit. Behavior of 
the evolutions of the couplings at large t reflects the effects of the naive di¬ 
mensions di ^ of the couplings which appear the naive dimensional part (|l^. 
These couplings evolute approximately proportional to exp{di^t) at large t, 
e.g., (72,0 ~ exp (2t), ( 74,0 ~const., (76 ,0 ~ exp(—2t). However, we remark that, 
as is seen in Eig. g, the couplings with a negative dimension may increase in 
the small t region, by the effects of the complicated interactions included in the 
bosonic, the fermionic and BE parts in the r. h. s. of Eq. ([T^). This means that 
these couplings may be important if we interested in the ” high energy” physics 
whose energy scale is close to Aq. |, |T1,||T^|T) 

As is seen in §2, the BE part does not contribute directly to the evolution 
of the coupling Ci^q. In Figs. || and we show the contributions of the bosonic 
and the fermionic parts in the r. h. s. of Eq. ( p^ for the evolution of ( 74,0 
coupling. The initial conditions are the same as in Figs. The bosonic 

part and fermionic parts cancel each other only partially, and they yield some 
amount of the total contribution. (See Fig. ^.) In Figs. |^|ll|, we show the 
contributions of the bosonic, the fermionic and the BE parts for the evolution of 
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Ci,i coupling. The initial conditions are the same as above. From the figures, 
we see that the contributions of the bosonic and the fermionic parts hardly 
cancel the dominant BF part. However, the whole contribution is much smaller 
than in the case of € 4 ,^. (See Fig. p^ .) 

To see the approximate discrete 75 -symmetry, in Fig. |^, we show the abso¬ 
lute value of the couplings Ci^ at t = 4 where almost all quantum corrections 
have been included into the effective generalized potential. In this figure, we 
have dividedized the coupling by the factor 




( Qexp (t) \ i 

V Ao ; 


(26) 


where gs = —Ci^i and Q is a typical energy scale of the system, due to the 
following two reasons. The first, the coupling C/,o is approximately propor¬ 
tional to the factor exp (d^pt) = exp {(4 — l)t} as is seen above. Therefore, we 

have divided the couplings by which is proportional to exp{(4 —Z)t} 

to eliminate the exponential t-dependence of the couplings. The second, the 
evolution of the coupling Ci^o is the order C{ i at each step of the the NPRG 
equation, since the equation has the one-loop structure. For example, Cs^o is 
the order of (Ci^i)^ as is seen in Fig. E Therefore, we have divided the cou¬ 
plings by \C\i \ = g\. Note, this normalization is equivalent to normalized the 
original couplings without the caret by the factor Q‘^~^g\. In calculations, we 
put Q/Aq = 0.5. In Fig. |^, we see that the odd couplings are somewhat sup¬ 
pressed in comparison with the even couplings in the neighorhood, except for 
the coupling Ci^o which is not small. 

It seem that approximate 75 -symmetry is hold to some extent through the 
evolution except for Ci^o- However, nonvanishing value of the coupling Ci^o 
induces the important effects. In fig. |^, we show the generalized potential U 
at CT = 0 when t becomes 4. It is seen that U does not have minimum at </) = 0 
because of the nonvanishing Ci^o- The cr-meson field (j) has non-zero vacuum 
expectation value < >. Therefore, we must shift the 0-field to the minimum 

to consider the real particle mass of the a-meson. This shift may enlarge the 
other odd couplings as in the case of the spontaneous symmetry breaking (SSB) 
and make the nucleon mass larger, since Ci^i < (f) >= —gs < 4 > > is positive. 
However, we remark that this enlargement of the odd couplings Ci^ {I > 3) is 
caused by the nonvanishing Ci^o and not by the negative C 2 ,o as in the case of 
SSB. This enhancement of the symmetry breaking may yield the possibility for 
the following scenario. 

(1) In the initial bare action, there exists approximate 75 -symmetry and the 
nucleon mass is not so large at t = 0. 

(2) Quantum corrections yield the large nonvanishing Ci^o and the field (j) 
has a non-zero vacuum expectation value which causes the redefinition of ip. 
This shift of field (p enlarge the other odd couplings and yields the large mass 
of nucleon. 
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In §5, we examine this scenario phenomenologically using the Hartree ap¬ 
proximation. 


4 Relation between the local-potential 

approximation and the Hartree approxima¬ 
tion 


In this section, we discuss the relations between the local-potential approxi¬ 
mation in the NPRG equations and the traditional Hartree approximation in 
quantum hadrodynamics. 0110^0 We truncate the r. h. s. of Eq. 

according to the observations obtained by solving the NPRG equations. 
First, we perform the inverse replacement of Eqs. and (|^) (see also (p^), 
i.e., perform the replacement 

kr=V^, H = (l> A3d = (T = Q,™ (27) 


where we restored the caret notation defined in the beginning of the previous 
section. By these replacements, the naive dimensional part in Eq. ( |l^ vanishes 
and we get 

where the effective potential V is defined by 


l,m—0 


and 


n{(j), cr; A) 


r2((^, <7]t) = 2 


aV„ 

p2 + y2 





(29) 


(30) 


We remark that the fields and the couplings without the caret in these equations 
are the ones after the replacements ( ^ ) and have the ordinary dimensions. The 
f{A,n) in Eq. (||) is given by 

/(A, /i) = In (/iVA^) - ^ In (mVA^) (31) 

where p is an arbitrary scale parameter with the dimension of mass. The /(A, p) 
does not depend on the field (j) and a. Therefore, this constant term contributes 
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only to the evolution of the coupling Cq.o which is the absolute value of the 
generalized potential at the origin (ti, cr) = (0,0). Since Co.o rnust be deter¬ 
mined phenomenologically in the actual calculations, we drop /(A,/r) below for 
simplicity of the notation. 

Next, we divide V in three parts. 

F = (</>; A) + (</>; A)a + VMgher .(<(-, a; A), (32) 

where 


y(°)(<^;A) = yW(<^; A) !<(>', 

1^0 1^0 

and VMgherA(l>, a; A) = ^ ^ Ci,^4>^a^. (33) 

m—2 l—O 

In Fig. 1^, we show the tr-dependence of the effective potential V at i^/Aq = 
0.1 and t = 4. The initial conditions are the same as in the previous section. 
It is seen that V is well approximated by the linear function of a. Therefore, 
below, we neglect the higher a terms Uhigher a- By this approximation, Eq. 
(p^ is reduced to following two equations. 


-a-'^e(°)(<(>;A) 


-a-^^fW(</>;A) 


ln{(p2 + V'i°^(</,;A))//i2} 


IOtt^ 


— In{(p2 + (FW(0;A))W} 


(34) 


1 


167r2 


<^(</';A) 


p 2 + E(°)(<^;A) 


^I/«(^.A)(eW(0,A))2 

1 

{P^ + <^(<(>;A)}{p2 + (EW(0;A))2}' 


(35) 


We remark that the BF part of Eq. (H) dropped in the Eq. (|34|) as is seen in 
the §2. 

Equations ( p4|) and (^) are coupled differential equations for the potential 
A) and A). However, these differential equations are difficult to 

be integrated out since the unknown A-dependence of the couplings appears 
in the r. h. s. of the equations. Therefore, we divide the integrations into 
two steps. First, using the original renormalization group equation (IH), we 
calculate the coupling evolutions from t = 0 to t = t* = In (Aq/A*) below which 
the couplings Ci^m vary slowly. In figs. 0 and |T^ we show the t-dependence 
of the couplings C' 2 ,o/Aq and Ce.oAg after the replacement of (p^). In the large 
region of t, in fact, they vary slowly. Therefore, below t*, we perform the 
integration by putting the couplings in the r. h. s. in the differential equations 
to be constants. IT^ 
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For the constant couplings, we use the couplings at A = 0 rather than ones 
at A = A*. This approximation may be justified by the following reason. For 
example, consider the integral 


r dAA3ln{(/ + yW(0;A))^2}. 

Jo 


(36) 


We expand this integral in powers of (j) around (p = 0. 




n! i9(/>” 

/ dAA^\n{{p‘^+ vll\(l)]A))/fj.‘^}\^=o 

Jo 


n—0 
rA* 


+ E 


1 a" 




—^ n! dp' 


dAA-^ 


p2 + y(°)(0;A) 




(37) 


In actual calculations, we must determine some lower terms of p phenomeno¬ 
logically. For example, in the ordinary renormalization procedure, in we must 
determine the coefiicient of when n is smaller than 5. In the case of the finite 
cutoff field theory, we may need to determine some higher terms to remove the 
cutoff-dependence of the physical results. Therefore, the quantities 

we can calculate are the coefficients of the higher terms in Eq. (|^). In such 
a term, the integrand has large power of p in the denominator. This means 
that such a integral is dominated by the low-energy behavior of the integrand. 
Therefore, we use the low-energy couplings at A = 0 rather than the high-energy 
couplings at A = A*. 

Using the couplings at A = 0, we get 


I/(°)(</);0) = U(0)(</,;A*) + ^^^ dAA3ln{(p2 + ujy(</);0))//r2} 

- dAA3ln{(p2 + (U«(0;O))2)//r2} 

= V^^\p-A*) + \ J ^0(A*-b|)ln{(p^ + U®(0;O))//r2} 

- IJ - bl)Tr[4 ® J,]ln{(p^ + (UW(<(>;0))2 )/m^} 

(38) 

U«(</>;0) = V^^\P;A*) + ^,vil\p-,0) r dAA^ J 

iOTT Jo p2 + y w 0) 

- £ dAA^ 
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X 


1 


{X + v;7(</';0))}{p^ + {^w(</';0))^} 

- \p\)tA{p) 


(-l)dV 


i2ny 




i27ry 


0(A* - Ipl)iA(p) 


TrliGjp)] 
4A ’ 


where 


and 


A{p) = 


-1 


p^ + vil\y-,o) 

■ p0)ls 


G{p) = 




(39) 

(40) 

(41) 


p2 + (T/(l)(<^; 0))2 

are the cr-meson and the nucleon full propagators, respectively, Ig is an unit 
matrix in spinor space, and the traces are taken over both the spinor and isospin 
indices. 

The first terms in the r. h. s. of Eqs. (|3|) and (^) play a role of counter¬ 
terms in the ordinary renormalization procedure. Second and third terms of the 
r. h. s. of Eq. ( p^ are the bosonic and the fermionic one-loop contributions, 
respectively. Equation (39) is an equation for calculating a scalar part of the 
nucleon self-energy, i.e., an effective nucleon mass. Especially, the third term in 
the r. h. s. corresponds to the Fock-term. In Figs, [l^a) and (b), we show the 
graphical expressions of each terms in Eqs. (^) and (|39[). 

If we consider the second derivative of the both sides of Eq. (^) with respect 
to (j), we get the following equation for calculating the effective cr-meson mass 
which is defined at the point where external momentum vanish. 


Vll\y;0) = K7(0;A*) 


( 0 ), 


4>4> 


1 


1 

+ 2 


{-l)d^p 


i2ny 

{-l)d*p 


- bl)*A(p) 


J - \p\mip)r 


(-zE«(</.;0)) J - b|)Tr{*G(p)}(-l) 

{-l)d‘^p 




(271-) 


4 0(A* - b|)Tr{*G(p)^G(p)}(-l) 


(42) 


The graphical expression of each term in Eq. (|42|) is shown in Fig. |20|(a). 
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Similarly, If we differentiate the both sides of Eq. with respect to (j), 
we get the following equation for calculating effective Yukawa coupling with 
vanishing external momentum. 




X 


X 

+ 


A*) + J - \p\Wp) 

2 (-rF«(</.,0))HY«(<^;0)) J 


- bl)*A(p) 


Tr{zG(p)} 

4A 


(-.V'y>w,0))»(-.r^« W;0)) J 

m- - H)(.A(rt)^Al!® 

(-.ry>W;0))» / 


4A 


(43) 


The graphical expression of each term in Eq. ( ^ ) is shown in Fig. ©b). Using 
the Eq. (B^, we can calculate the effective Yukawa couplings. 

In fig. |21|, we show the (/)-dependence of which is obtained at t = 4 by 
solving the NPRG equations. The initial conditions are the same as above. We 
see that can be well approximated by a linear approximation 


(0; A) = Co,i + = M' - (44) 


where M' = Cop. 

If we put Eq. (^) into Eqs. ( ^ and regard V^^{4>]Q) as a constant pa¬ 
rameter in the r. h. s. in Eq. (|^), we get the usual Hartree approximation.^, 
II, 10, since the bosonic loop in Eq. (^) yields only the constant contri¬ 
bution which is determined phenomenologically and it can be dropped. Under 
the same approximation, the second term of r. h. s. of Eq. ( ^^ is dropped, 
Eq. (|3^) reduces to the equation for calculating the Fock-term contribution 
with vanishing external momentum, using the Hartree propagator pertur- 
batively. Similarly, only the first and last terms of r. h. s. of Eqs. (^2|) 
and (^ ) remain, and we get the equation in the usual RPA|^, ^ and vertex 
correction||llj, respectively. 

The results presented here is the ones at zero baryon density. To justify 
the Hartree approximation at finite baryon density, it is needed to extend the 
NPRG equations to finite density. Furthermore, the effects of the other mesons, 
especially of the w-meson should be examined to study the nuclear matter. 
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However, below, we assume the validity of the Hartree approximation at finite 
density in the a-Lo model. In fact, according to the calculations based on the 
Nuclear Schwinger-Dyson formalism, [Q ^ 

corrections beyond the Hartree approximation seem to cancel each other. [p0[ In 
§5, we examine the approximate 75 -symmetry in the Hartree calculations in the 
a-uj model at finite baryon density. 


28 


the higher order 


5 Approximate 75 -symmetry in the Hartree 
calculations in nuclear matter 


In this section, we examine an approximate 75 -symmetry in u-lo model 1^, ^ in 

the framework of the finite cutoff field theory, 0 , 1 , 1 , 0 ,|T| using the Hartree 
approximation. In this model, the energy density of nuclear matter is 

given by 

e(p,$,A0 = t^($) + Ka™(d>;A0-zy'^Tr[Gz5(p)] + ^p2. 

(45) 

1^0 

A') = -I j ® h] ln{(p" + (47) 

In these equations. A' is a convenient cutoff which is introduced by hand for 
calculations and corresponds to A* in the previous section. M is a nucleon mass, 
and 4), c;, M* and Cy are defined by 

$ = 5 < 0 >, Ci = Cio M* = M-^, and Cy = (48) 

’ iriy 


, respectively, where niy and gy are w-meson mass and w-nucleon coupling, 
respectively, and < (j) > is the ground state expectation value of (j) at each 
baryon density. The baryon density p is given by 


P = 



(49) 


where is a Fermi momentum. Gd{p) is the density part[0 0,0] of the nucleon 
propagator which depends pp explicitly, and is given by 


Gd{p) = {-i^-p + M*h)27r6{p*^+M*^)e{po)e{pF-\p\)<»Ii; 

Pp = (7'4 + *VFo,p); Po=*P4, (50) 
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where p are the three dimensional momentum and Wo is given by 


Wo = gv < Wo >= (51) 

with the ground-state expectation value of the time component wo of the co- 
meson field. $ is determined by the equation of motion of the cr-meson, i.e., 

^e(p,$;A')=0. (52) 

The first and the second terms in the r. h. s. of ( ^5| ) correspond to the first 
and the third terms in the r. h. s. of (^). The first term plays a role of the 
’’counter-term” in the ordinary renormalization procedure. We have truncated 
the higher-terms {I > 6 ) in it. Coupling c; is determined by the condition 


1 . n 4 


(53) 


where c* are effective or ” physical” couplings which are determined phenomeno¬ 
logically. We choose 


c 

9 


* 

I 


S 


0 (/ = 0,1,3,4), 


and 


C2 


(54) 


where Cs = gsM/iris and rus is a cr-meson mass. Note that the model depends 
on iris, mv, gs, gv and c/ only through the raitos gs/rris, gv/rus and cifgs- If 
we use the nucleon mass M = 939MeV, only four variable parameters remains, 
i.e., Cs, Cd, co/gs and Colg\- We choose the four parameters to satisfy the 
saturation and compressional properties in nuclear matter. |p^ For the satura¬ 
tion properties, we use the binding energy 15.75MeV at the saturation density 
Pst — 0.15fm“^. For the compressional properties, we use the results analized 
by Pearson, In table I, some available parameter sets are shown. We remark 
that the phenomenological determinations of C 5 and cq make the A'-dependence 
of the physical results very small. |]^ 

First, for each parameter set, we expand the energy density in powers of $ 
around $ = 0 ( or M* = M )at zero baryon density as follows. 

e{p = 0, $; A') = ^ i^e(p = 0, $; A*)W=o<^^ = E (55) 

In Figs. Hjand we show the absolute values of the couplings /g'-^ (1 < 

I < 6 ) which are determined phenomenologically in this model as is seen above. 
We use the nucleon mass M as the typical energy mass scale Q which is used 
to normalized the couplings. In each case of the parameter set (PS), it seems 
that there are no explicit suppression of the odd couplings. 
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Next, we expand the energy density around M* = M*^ at the saturation 
density pst, i.e., 


1 

£{p = A') = ^ = Pst,^]h')\<s>=<f>, 

1=0 


‘j’' = E 

1=0 


(56) 


where $st = M—M*^ is the expection value of the cr-meson field at the saturation 
density psf The absolute values of these couplings are shown in Figs. ^ and 
p5| . It seems that there is an suppression of the odd couplings b*. There may 
be the approximate discrete 75 -symmetry at finite baryon density. 

Finally, we expand e(p, $; A') around M* = M*^ but put p = 0, i.e.. 


s{p = 0 , $;A') = 


1=0 




1=0 


gi 


(57) 


The absolute values of these couplings are shown in Figs. and We see 
that the coupling O;* /g^ have similar features as in the case of p = pst, except 
for 1 = 1. Naturally, nonvanishing a* corresponds to the fact that $ = $st is 
not the solution of the equation of motion (52) at zero density. 

Although we did not show the results, the parameter set PS9 in Ref. 10) 
which has also small M*^ (= 0.5116M) at pst yields the similar results as in 
the cases of PSIO and PSll. On the other hand, the parameter set PS4 in Ref. 
10) which has large M*^ (= 0.8887M) at pst yields no suppression of the odd 
coupling. However, the couplings C 5 and cq in PS4 are too large in comparison 
with the ones in PSIO and PSll, and it may not be the realistic parameter set. 

From these results, it seems that, there may be a approximate symmetry 
even at p = 0 but nonvanishing a* causes the sift of the field cj), yielding the 
enlargement of the other odd couplings and the large nucleon mass M. At the 
saturation density, the density contribution cancels the large al and makes the 
nucleon mass smaller. 


6 Summary 

In this paper, we have studied the quantum hadrodynamics using the non- 
perturbative renormalization group equations. The results obtained here are 
summarized as follows. 

( 1 ) At zero density, the NPRG equations are formulated in the quantum 
hadrodynamics (cr-nucleon model). 

(2) Approximate discrete 75 -symmetry is studied. It is found that the ap¬ 
proximate symmetry is retained to some extent through the whole evolutions in 
analogy of the exact one. However, the evolution makes the linear (/)-term large 
which causes the sift of the field (j), yielding the enlargement of the other odd 
couplings and the large nucleon mass M. 
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(3) Relations between the local-potential approximation in NPRG equations 
and the Hartree approximation are studied. The importance of the one-loop 
contribution is shown. The local-potential approximation contains not only the 
Hartree contribution but also the contributions of the Hartree-Fock, the RPA 
and the vertex corrections with vanishing external momentum. 

(4) The approximate discrete 75 -symmetry is studied phenomenologically 
using the Hartree approximation. There may be the restoration of this sym¬ 
metry at finite density, since the medium effects cancels the large linear (jj-term 
and makes the nucleon mass smaller. 

There are some interesting extensions of this work. 

(1) To formulate the NPRG equations at hnite baryon density/finite tem¬ 
perature. 

(2) To include the w-meson to study the saturation mechanism. 

(3) To include the 7 r-meson to study exact and approximate chiral symmetry. 

They are now under the studies. 
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K 

M*JM 


c'i 

CbM/g^, 

CQM'^/g'l 

PSIO 

300.0 

0.5156 

378.50 

288.38 

0.0024938 

-0.0018217 

PSll 

350.0 

0.5442 

358.11 

270.07 

0.0025532 

-0.0018579 


Table I; Parameter sets. is the effective nucleon mass at the saturation 

density. K is the incompressibility of nuclear matter and shown in MeV. To get 
these parameter sets, we have used A' =1.5GeV. 
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Figure 6 : Bosonic contribution for the evolution of the coupling (74^0 
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Figure 8: Total interaction contribution for the evolution of the coupling 
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Figure 10: Fermionic contribution for the evolution of the coupling Ci^i 
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Figure 14: The one-loop contribution with the Yukawa coupling Ci^i for 
The solid and dotted lines are fermion and boson lines, respectively. 


34 





Figure 15: Effective generalized potential at ct = 0 and t = A. 
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Figure 16: V at (p = O.IAq and 1 = 4. 
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(b) 


Figure 19: The graphical expression for ((a))and ((b)). 
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Figure 20: The graphical expression for ((a)) and ((b)). 
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Figure 21: at t = 4. 
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Figure 23: The absolute value of c*M* ^!g\ in PSll. 
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